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^ : ABSTRACT 

! Fast X-ray transients (FXTs) with timescales from seconds to hours have 

^ ! been seen by numerous space instruments. Because they occur at unpredictable 

! locations, they are difficult to observe with narrow-field instruments. Only a few 

^ ! hundred have been detected, although their all-sky rate is in the tens of thou- 

! sands per year. We have assembled archival data from Ariel-5, HEAO-1 (A-1 

^ '. and A-2), WATCH, ROSAT, and Einstein to produce a global fluence- frequency 

^ ! relationship for these events. Fitting the log N-log S distribution over several or- 

'. ders of magnitude to simple power law we find a slope of — 1.0^9 3. The sources of 

O ■ 

o ■ 
^: 

i their individual luminosity distributions are both fiatter than the -3/2 power law 

j> I that applies to uniformly distributed standard candles. Extrapolating from the 

^ I BATSE catalog of GRBs, we find that the fraction of X-ray fiashes that can be 

^ I the X-ray counterparts of gamma-ray bursts is a function of fluence. The exact 

fraction of GRB-induced X-ray counterparts is sensitive to the Rxf-y distribution, 
which we estimate from available experimental measurements. Certainly most 
FXTs are not counterparts of standard gamma-ray bursts. The fraction of FXTs 
from non-GRB sources, such as magnetic stars, is greatest for the faintest FXTs. 
Our understanding of the FXT phenomenon remains limited and would greatly 
benefit from a large, homogeneous data set, which requires a wide- field, sensitive 
instrument. 

Subject headings: gamma rays: bursts — stars: activity — stars: late-type — 
X-rays: binaries — X-rays: bursts — X-rays: general — X-rays: stars 



FXTs are undoubtedly heterogeneous, representing several physical phenomena; 
the a ~ -1 power law is an approximate result of the summation of these multiple 
sources. Two major contributions come from gamma-ray bursts and stellar flares. 
These two types of progenitors are distributed isotropically in the sky, however 
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1. INTRODUCTION 

Essentially all missions sensitive to cosmic X-rays have detected intense X-ray out- 
bursts with timescales from seconds to hours. Peak fluxes range to 1 Crab and above 
(1 Crab=2xl0^^ ergs cm^^ s^^ in 2-10 keV band). X-ray outbursts are unlike classical X- 
ray transients, which persist for weeks or months. X-ray outbursts fade within a day and are 
normally only seen once. Historically, X-ray outbursts seen for less than a day that lack per- 
sistent X-ray counterparts were called "fast X-ray transients" (FXTs). Recently a new term, 
"X-ray flashes" (XRFs), was coined for short (less than 1000 sec) intense bursts of X-ray 
emission (Heise et al. 2001). XRFs lack detectable gamma-ray emission, but are reminiscent 
of gamma-ray bursts in their time history. It can be diflicult in some cases to distinguish 
XRFs and other FXTs using archival data e.g. where the time resolution is inadequate. 

The sources of FXTs might be heterogenous; the experiments that measure them are 
certainly so. FXTs have been surveyed in a variety of energy bands, sometimes with low 
spectral and time resolution, and typically with low angular accuracy. There is no straight- 
forward way to obtain a uniform large set of FXT parameters. It is no surprise that diflerent 
experiments have suggested different sources for FXT, including flare stars, compact objects, 
extragalactic sources, and X-ray emission from GRBs [e.g. Pye & McHardy 1983 (hereafter 
PM83); Ambruster & Wood 1986 (hereafter AW86); Castro-Tirado et al. 1999 (hereafter 
CT99)]. However, the quantitative contribution of these several sources is still unknown, 
and room remains for unknown phenomena. The discovery by BeppoSAX of long-lasting 
X-ray emission from GRBs, the so called X-ray afterglow (Piro et al. 1998), and the strong 
prompt X-ray emission from GRB seen by Ginga (Strohmayer et al. 1998) and WATGH 
(Sazonov et al. 1998), shows that some fast X-ray transients are related to GRBs. 

Theoretical work predicts the existence of GRB types that are X-ray rich or gamma- 
ray quiet. This might be caused by physical properties of the GRB central engine and 
surrounding media, by the geometry of GRB source with respect to the observer's line of 
sight, or by the redshift of classical (standard) GRBs (see recent review by Meszaros 2002). 

In this paper we attempt to construct a uniform data set of FXTs from the heterogeneous 
data available to date. In section 2 we describe the experiments, and in section 3 our data 
analysis. In section 4 we discuss the connection of FXTs and GRBs, and estimate the fraction 
of GRB-related FXTs in FXT data. In section 5 we discuss other possible sources of FXT 
progenitors. 
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2. 



EXPERIMENTAL DATA 



The importance of sensitive uniform sampling has been well recognized in the case of 
GRBs (see Lee & Petrosian 1996; Petrosian & Lee 1996, and elsewhere) . The BATSE/CGRO 
experiment, with its large homogenous and uniform database, has allowed the study of sta- 
tistical properties of classical GRB. A similar database, if it existed, would allow an unbiased 
determination of the important parameters of FXTs, and could test different hypotheses for 
their origin. Although FXTs are by no means rare, we still do not have a large uniform set 
of these events. In its absence, we have combined the data provided by several very different 
experiments to quantify, as best we can, the statistics of FXT occurrence. While we recog- 
nize the imperfection of this approach, it is the only way to estimate the FXT distribution 
over several orders of magnitude in fluence. A brief description of the instruments and data 
sets is given below. 



PM83 published a catalog of FXTs detected in the 2-18 keV band by Ariel-5. Flux 
measurements from the rapid-spinning satellite were integrated in 100 minute bins corre- 
sponding to the satellite orbit. These fluxes provide a good estimate of the integral fluence 
in the time bin. Transient events were selected by their detection in a single, or at most a 
few, sateUite orbits. The hghtcurves therefore have 100-min time resolution. The limiting 
sensitivity was ~20 mCrab (5.5(t) in 100 minutes. PM83 found a cumulative log N-log F 
peak flux relation for 27 detected FXTs consistent with a power law with index a — -0.8 
± 0.5. PM83 estimated the total number of FXTs above their threshold as N=150-180 per 
year for the whole sky. 



The HEAO-1 satellite carried two X-ray surveys, Al and A2. These instruments scanned 
a great circle on the sky every 35 minutes. Sources were sampled for a short interval once 
or twice per scan (10 sec for Al and 60 sec for the A2 experiment). 

AW86 published a survey of transients detected by HEAO-1 Al in the 0.5-20 keV band. 
This AW86 survey was based on excess emission over an 12-hour interval. 10 FXTs were 
detected, with a hmiting sensitivity about ~4 mCrab. The log N-log F distribution was 
consistent with a power law with index a = — 1.0± —0.5. AW86 estimated the total number 
of FXTs as N=1500-3000 per year above their threshold, assuming that a typical event lasted 
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Ariel-5 
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more than 1.5 hours. 

Connors et al. (1986, hereafter C86) pubhshed the HEAO-1 A2 survey of FXTs in the 
2-20 keV band. 8 FXTs were detected with a hmiting sensitivity about ~4-6 mCrab. The 
log N-log F distribution fit a power law with index q;=-1.0±0.7. The duration of the detected 
FXTs was poorly determined, and might have anywhere between 60 s and ~2000 s. The 
all-sky rate of such FXTs would have been between 10^ per year for durations 2000 s and 
3x10^ per year for durations of 60 s. 

2.3. WATCH 

CT99 published the detection of 7 FXTs by the WATCH all-sky monitor onboard the 
Granat sateUite. FXTs were detected in a 8-15 keV energy channel. WATCH, with a time 

resolution of several seconds, was able to resolve the fast transients it detected (~100 min 
for 3 FXTs and '^1 day for the rest). The typical peak flux for these sources was several 
hundred mCrab. Because the lightcurves are well-sampled, we have a good estimate of the 
total fiuences. Sazonov et al. (1998) used WATCH to observe an additional 95 events in 
the energy band from 8 to 60 keV. 86 of them were observed by other GRB experiments, 
and confirmed to be GRBs, and 47 events were locahzed to 1 degree or better. We used the 
Sazonov et al. (1998) data to constrain the Rxj-y ratio for gamma-ray bursts (see section 4). 

2.4. ROSAT 

Vikhlinin (1998, hereafter V98) searched for faint X-ray bursts (duration 10-300 s) in 
pointed ROSAT PSPC observations with a total exposure of 1.6x10^ s. Only the 'hard' 
ROSAT energy channel (0.5-2 keV) was used for their study. A total of 141 X-ray fiares 
with duration 100-300 s were detected. The quiescent fiux was consistent with zero in ap- 
proximately half of the burst sources, while positive quiescent emission was detected for the 
rest. 112 bursts were identified as Galactic stars. One burst was from the X-ray binary LMC 
X-4. The 28 remaining FXTs were not identified, but for most of them starlike optical coun- 
terparts were found in the Digital Sky Survey. The limiting fiuence of these events (10-100 s 
duration) corresponds to 2.6x10^^° ergs cm~^ in the 0.5-2 keV band. We have translated 
the fiuence from ROSAT energy band to 2-10 keV under assumption of Crab-like spectrum 
and neglecting foreground absorption. This implies a limiting fiuence of 3.5 mCrabxsec. 
The estimated rate of such events is ~10^ per year for the whole sky. 

Greiner et al. (2000) searched for GRB afterglows in the ROSAT all-sky survey. During 
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the survey, the telescope's field of view scanned a great circle on the sky, sampling sources for 
10 to 30 s. The sensitivity limit per scan is 10~^^ erg s~^ cm~^. Greiner et al. (2000) applied 
tight restrictions on the time profile to extract GRB afterglow candidates, and selected 23 
afterglow candidates. Closer examination showed that about 14 of afterglow candidates were 
in fact fiares from late-type stars of class M-dMe. With less restrictive requirements on the 
time profile Greiner et al. (2000) found more afterglow candidates, mostly associated with 
flare stars. 



2.5. Einstein 

Gotthelf, Hamilton & Helfand (1996; hereafter GHH96) reported the detection of 42 
faint X-ray flares in the Einstein IPC data. Those flares have typical fluxes from 10"^*^ ergs cm" 
to 10"^ ergs cm^^ in the 0.2-3.5 keV energy range, and their durations were less than 10 s. 
Their detection rate corresponds to about 10^ flares per year for whole sky. However, they 
concluded that only the hardest of their events, consisting of 18 flashes, was undoubtedly of 
cosmic origin. The remaining soft events might have been detector artifacts. 

2.6. Ginga 

Strohmayer et al. (1998) reported on spectral flts and the Rxf-y ratio for 22 GRBs 
observed with Ginga over 2 to 400 keV band. We used their flts, together with data from 
other instruments, in our Rx/j distribution. 

2.7. BATSE 

BATSE GRBs are cataloged in the 4B and 'Current' BATSE databases^. A recent 
catalog of bright BATSE GRBs was pubhshed by Preece et al. (2000, hereafter P2000). 
Two searches by Kommers et al. (1999) and Stern et al. (2001) showed that the BATSE 
database contains almost as many untriggered GRBs as triggered GRBs. The untriggered 
GRBs tend to be fainter than the triggered GRB (see Fig.8 of Stern et al. 2001), for natural 
reasons. We used both triggered and untriggered GRBs to obtain a log N-log S for GRBs in 
the gamma-ray band. 



^http://gamniaray.nisfc. nasa.gov/batse/grb/catalog/ 
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Spectral fits for 53 GRBs detected with BATSE were published by Band et al. (1993). 
Using the P2000 catalog, we calculated spectral parameters and estimated Rx/i for 81 more 
BATSE GRBs. 



2.8. BeppoSAX 

We used BeppoSAX afterglow observations (see review in Prontera et al. 2000a) as in- 
puts to for our estimate of the fraction of GRB-related FXTs. Prontera et al. (2000a, 2000b, 
2001), in't Zand et al. (1999, 2000a, 2001), Pian et al. (2001) and Nicastro et al. (2001) 
presented the X-ray and gamma-ray fluences for 15 GRBs in the 2-10 keV and 40-700 keV 
bands measured by BeppoSAX. We use these data to help estimate Rxj-y- Although we 
discuss the statistics of GRBs and XRPs reported by Heise, in't Zand, & Kuulkers (2000), 
and (Heise & in't Zand 2001; Heise et al. 2001) we did not include them in our quantitative 
analysis because no numerical data for those XRPs have been pubhshed in useful form. 

2.9. RXTE ASM and Konus/Wind 

Smith et al. (2001) recently presented RXTE ASM and Konus/Wind fluence measure- 
ments in the 1.5-12 keV and 50-200 kcV energy bands for 14 GRBs detected by both these 
instruments. We have used these data for our -Rx/7 estimate. 

3. DATA ANALYSIS 

Our study started from the analyses published for the individual instruments. We used 
both published results, and, in some cases, data from public archives. However, because 
results of previous studies were strongly instrument dependent, they cannot be directly 
compared. We need to convert data from the individual instruments into a quasi-uniform 
data set. 

Our analysis will proceed in two directions: in the forward direction (approach 7^1), 
we will fit a power law luminosity distribution (log N-log S) to the data. In the backward 
direction (approach #2) we extract a model-independent luminosity distribution by plotting 
each data set on the log N-log S plane. Approach #2 requires the several assumptions 
discussed below. 

Most of the aforementioned experiments measured the fluence (time-integrated flux) of 
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detected events, but the HEAO-1 instruments provided only intermittent sampling of fluxes. 
These HEAO-1 results are important, however, because they were able to detect fainter 
FXTs than Ariel-5 and WATCH. ROSAT and Einstein detected even fainter FXTs, but 
those experiments were sensitive in a much softer energy band, and are subject to additional 
corrections when converted to the 2-10 kcV band. The WATCH data also need correction, 
because that instrument was not sensitive to softer X-rays, but in this case the overlap with 
other X-ray instruments is larger and hence the correction is not so critical. Additionally, 
interstellar absorption, which can be safely neglected in most cases above 2 keV, is an 
important factor in the ROSAT and Einstein bands, adding more uncertainty to the results. 
Finally, we did not have information about individual FXTs detected with ROSAT and 
Einstein, and used instead the integral estimates presented by V98 and GHH96. 

For our primary analysis (approach #1 above) we have used data from Ariel-5 (en- 
ergy range 2-18 keV, minimum detectable fluence about 100 Crab x sec), HEAO-1 Al & 
A2 experiments energy range 0.5 - 20 keV & 2 - 20 kcV, minimum detectable fluence 0.5 
Crab X sec), and WATCH (energy range 8-15 kcV, minimum detectable fluence about sev- 
eral hundred Crab x sec). The fluence sensitivity limits above are given for long-duration 
events, up to 1 day. The fluence sensitivity is much higher for short intense bursts, where 
background is not important. 

To include HEAO-1 data into our analysis, we need to make assumptions about the 
duration and time profile of the transients. Fortunately, because of the nearly ~1 slope 
of the cumulative distribution, almost the same distribution follows from a wide range of 
assumptions on duration. We have only lower (the scan duration) and upper (the scan 
interval) limits to the duration of the HEAO events. We estimated fluences for the HEAO-1 
events by multiplying the measured flux by an assumed event duration. Our initial estimate 
of event duration was taken as the geometric mean of the minimum possible FXT duration 
(usually the scan duration of 10-60 sec) and the maximum possible FXT duration, which is 
several hours for the Al experiment and ~4000 s for the A2 experiment. 

To quantify the transient statistics, we assume that they are distributed as A^(> S) — 
Aq X 5""", where Ao is the normalization at 1 Crab x sec, and determine the parameters of 
this distribution with a maximum likelihood analysis. In our first analysis, we assumed that 
the spectrum of the events was similar to the Crab nebula, i.e., a power law with photon 
index F = —2. We then calculated the minimum detection fiuence for each experiment. The 
overall event rates were taken from the published results of the original analyses (AW86; 
C86; PM83; CT99). 

More detailed discussion of sky coverage, thresholds and detection bands can be found 
in the references above. The parameters of our maximum-likelihood determination of the 
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parameters a and are shown in Fig.l. 

It is known that many flare stars produce short, intense X-ray flares (Strassmeier et al. 1993) 
with thermal spectra. Because many FXT have been identified with fiare stars (AW86; C86), 
we have calculated the FXT distribution that results if the FXT consist of mixture of 
Crab-hke and thermal spectra spectrum. A mixture of 30% thermal-spectrum FXT, with 
kT — 0.7A;ey, and 70% Crab-spectrum FXTs maximizes the maximum likelihood function, 
considering the range kT = 0.5 - 4.5 keV, and 'thermal' burst fraction to 100%. 

Given this mix of FXT spectra, our best estimate of a^A^ yields q;=0.95±0.1 and 
= (4 ± 1) X 10^ at ICrab x sec. A conservative, 3a limit bounds the parameters to 
q;=[0.7-1.2], Ao=[10^-10^] In any case, FXTs are common, occurring somewhere on the sky 
dozens to hundreds of times each day. They are one or two orders of magnitude more 
common than classical gamma-ray bursts (Kommers et al. 1999). Note that it is impossible 
to determine whether all the events that make up this distribution fit our strict definition 
of occurring only once. If some of the FXTs come from magnetic stars, RS CVns, or other 
nearby objects, they may indeed repeat. These repetitions would have been seen by HEAO-l, 
Ariel-5, etc. if those instruments had better sky coverage or longer lifetimes. 

Our distribution ( a = 0.95 and Aq = 4 x lO'', 70% Crab-like vs. 30% thermal spectra) 
is shown on Fig. 2. It is not significantly different from the distribution that follows if we 
assume that all bursts have a Crab-hke spectra. 

Fig. 2 also shows the segments of the distribution seen by the individual instruments (our 
analysis approach 7^2). The shaded areas represent la errors for the individual experiment 
log N-log S curves, including the normalization and the statistics of the number of events 
seen by each experiment. 

Although the duration of the HEAO-l events is quite uncertain, this uncertainty does 
not affect the overall distribution. We plot the backward-derived log N-log S assuming the 
geometric-mean time duration discussed above. But because the overall distribution has a 
slope close to 1, changes in the assumed duration simply shde the log N-log S contribution 
from the HEAO-l instruments up or down along the same curve. The dashed/dotted bound- 
aries show the effect of changing the duration to the extreme possibilities on the HEAO-l 
A2 log N-log S distribution. 

The ROSAT results from V98 and Einstein results from GHH96 were not used to deter- 
mine our net log N-log S distribution. However, they are plotted on Fig.2 for comparison. 

We have translated the fluence measured in the ROSAT energy band to 2-10 keV by assum- 
ing a Crab-like spectrum with negligible foreground absorption. The events reported by V98 
had fiuences from 3.5 to 200 mCrabxsec. We estimated the total number of events per year 
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for ROSAT events which had no detectable quiescent emission as ~2xl0^ year~^, and for 
events without optical counterparts as ~4xl0^ events year~^. We have taken these values 
as upper and lower limits to the FXT rate at this fiuence (shown as the cross on Fig. 2) The 
Einstein data were also converted to 2-10 keV assuming a Crab-like spectrum (shown as 
the rhomb on Fig. 2), which is consistent with the spectral shape found by GHH96 for hard 
events. 

Our integrated distribution is consistent with a power-law with a slope a = — l.Olo'g. 
A spatially homogenous distribution of identical sources (standard candles) would have a 
slope -1.5. Our distribution is consistent with the slope of ~0.8 derived by PM83 from 
Ariel-5 events. By adding data from several other experiments, we have extended the PM83 
distribution by several orders of magnitude towards the faint end. The ROSAT and Einstein 
points, derived from V98 and GHH96, are consistent with an extrapolation of our log N-log 
S curve. 

Fig.3 shows the sky distribution of the events used in Fig.2 (Ariel, HEAO, and WATCH 
marked with crosses, and ROSAT and Einstein events marked with triangles). The sky 
distribution is consistent with isotropy. The dipole moment of the distribution towards 
Galactic Center is 0.07±0.08 , and the quadrupole moment is 0.31±0.04 (for an isotropic 
distribution the dipole moment is zero, and the quadrupole moment is 1/3). We see no 
significant anisotropy for any subset of these data. However we need to note that there are 
no reliable exposure maps for mentioned above experiments, making it difficult to be definite 
about the lack of a disk population in these data. 



4. The Connection to GRB 

The connection between FXTs and GRBs is especially intriguing. As a result of the 
Ginga and WATCH studies of prompt X-ray emission from GRBs (Strohmaycr et al. 1998; 
Sazonov et al. 1998), and the BeppoSAX discovery of GRB afterglows (Piro ct al. 1998), it 
is evident that GRBs can emit a large fraction of their energy in the X-ray region (see 
also Frontera et al. 2000a). This emission is both coincident in time with the gamma- 
ray burst (prompt emission) and delayed (afterglow). Recently BeppoSAX has found a 
few GRB-like events with very little or no gamma-ray emission but strong X-ray emission 
(Heise & in't Zand 2001; Heise et al. 2001). These events were designated "X-ray flashes" 
(Heise & in't Zand 2001). However, Kippen et al. (2001, 2002) argue that they merge 
smoothly with regular gamma-ray bursts in a single distribution of their parameters. Such X- 
ray rich GRBs were estimated to total about 30% of all GRBs (Kippen et al. 2002). Whether 
XRFs and GRBs originate from a single class of progenitors with a smooth distribution of 
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parameters, or are at least in some cases intrinsically different, is not yet clear. Comparison 
of GRBs and XRFs is fraught with selection effects, including the fact that the observed 
XRFs were selected in a different energy band (X-rays) than gamma-ray bursts. 

In the absence of gamma-ray data, one might classify a conventional GRB that was 
detected in the X-ray band as an XRF. Grindlay (1999) examined Ariel-5 and HEAO-1 data 

for an excess of gamma-ray quiet FXTs, above those expected from conventional gamma-ray 
bursts, that might indicate beaming of GRBs. He estimated that the GRB share of Ariel-5 
and HEAO-1 FXT events was about 50%. We wish to refine this estimate by bringing several 
parameters into our analysis. These include the fraction of GRBs that emit in the X-rays, 
the ratio of the fluence emitted in the X-ray band to the fluence emitted in the gamma-ray 
band, the fraction of GRB that produce X-ray afterglows, and the ratio of the X-ray fluence 
emitted in prompt phase to the fluence emitted in the early afterglow (the first 10^ — 10^ 
sec) 

BeppoSAX and WATCH (Frontera et al. 2000a; Sazonov et al. 1998) indicate that the 
first parameter is close to unity for long-duration GRBs. Long-duration GRBs are defined 
as events with Tgo > 2 sec in the BATSE data set (Kouvehotou et al. 1993). The second 
parameter of interest is the ratio of the fluence emitted in the X-ray band (2-10 keV) to the 
fluence emitted in the gamma-ray band (50-300 keV), Rx/y The ratio Rx/-^/ can be estimated 
from the few events observed in both the X-ray and gamma-ray band, or, using the larger 
BATSE sample, by extrapolating the measured spectrum into the X-ray band. If the latter 
is consistent with the former estimate, one could use the full BATSE GRB database to find 
the Rx/-y distribution 

Band et al. (1993) has shown that GRB spectra can be fit by the expression: 

in the 20 keV-2 MeV energy band. Using typical values for the parameters (see for 
example Lloyd & Petrosian 1999; P2000), the average value of Rx/j is about several percent. 
However, Preece et al. (1996) has shown that such an approach leads to an underestimate 
of X-ray fluxes for about 15% of BATSE GRBs. The magnitude of this underestimate can 
be relatively high (up to order of magnitude above Band's formula) for an individual GRB. 

Ginga, which was able to simultaneously measure X- and gamma-ray emission from 
GRBs, found an average ratio of prompt emission Rx/j of 0.24 for 22 GRB events, although 
the logarithmic mean Rx/j was 0.07 (Strohmayer et al. 1998). Rx/j was close to unity for 
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several events. BeppoSAX confirmed that such a high ratio is not rare (Prontera et al. 2000a; 
Kippen et al. 2002). 

Data from different instruments are subject to significant differences in their energy 
bands and event selection criteria and often do not agree between themselves. To refine our 
estimate of Rxf-y we have used all relevant available data including BATSE, Ginga, RXTE, 
BeppoSAX and WATCH. Fits to the spectra of 53 BATSE and 22 Ginga GRBs, using the 
Band formula, are found in Band et. al (1993) and Strohmayer et al. (1998). Using the 
catalog of BATSE GRBs by P2000, we have calculated the Band parameters for 81 additional 
GRBs from BATSE data. For all these bursts we calculated Rx/-^ by integrating the Band 
formula. Smith et al. (2001) presented RXTE ASM and Konus measurements of 14 GRB 
in 1.5-12 keV and 50-200 keV energy bands. To convert ASM counts to X-ray fluence 
we fit the three-color ASM data to a power law for each GRB (assuming a single power 
law in the 1.5-12 keV region). To obtain the 50-300 keV fluence from the single-channel 
Konus data we assumed an average value for the Band parameters from P2000. BeppoSAX 
data for 15 GRBs are presented at Frontera et al. (2000a, 2000b, 2001), in't Zand et al. 
(1999,2000a,2000b,2001), Plan et al. (2001) and Nicastro et al. (2001), who give X-ray and 
gamma-ray fluences for the 2-10 keV and 40-700 keV bands. When available, we used the 
individual Band's parameters for each GRB, and otherwise used average values from P2000 
for 50-300 keV channel. 

Sazonov et al. (1998) published a catalog of 95 GRBs detected with WATCH. We used 
82 of these events, those for which fluences are available for both energy channels (8-20 keV 
and 20-60 keV). Full data were unavailable for the other 13 events. For the WATCH data 
we extrapolated Rx and R^ to our 2-10 keV and 50-300 keV channels by assuming that the 
a Band's parameter was fixed at average value found from Ginga, RXTE and BeppoSAX 
measurements, ~ -0.85, while the other parameters were fixed at average values from P2000. 

Fig. 4 shows the distribution of i?x/7 for all these events. The solid hne shows the 
total Rxi-y distribution derived from the Ginga, RXTE/Konus, BeppoSAX and WATCH 
measurements (134 GRB total), while the dashed hne represents the BATSE events (134 
GRB total). Note that the BATSE ratio comes from an extrapolation of the measured 
spectra to the X-ray band. The non-BATSE events, which are triggered at lower energies 
than BATSE, tend to have larger median Rx/-, - along with a tail of even higher Rxf-y, X-ray- 
rich events. We will call them X-ray detected GRB (XDGRB). The later distribution is in a 
good agreement with BeppoSAX measurements (Heise & in't Zand 2001). We need to note 
that the Rx/'y distribution certainly depends on how one weighs events from the different 
experiments. 

We consider the XDGRB estimate of Rxf-y to be more reliable than the 'BATSE' esti- 
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mate, because the 'BATSE' -Rx/7 ratio was calculated by extrapolation of Band's parameters 
from the 20 keV-2 McV region, while for the XDGRB Rx/-y was actually measured, or cal- 
culated using Band's parameters fit to 2 keV-700 keV measurements. The shift between 
the peaks of the two distributions is in the direction expected from selection effects (see 
Appendix) . 

To determine the integrated X-ray fluence from GRBs, we need to consider not only the 
prompt emission, but also the afterglow emission (the third and forth parameters mentioned 
at the beginning of this section). Data from BeppoSAX show that 80-100% of GRBs that 
emit prompt X-ray emission also show an X-ray afterglow (Frontera et al. 2000a). We use 
this fraction (100 %)as our estimate of the third parameter. 

To estimate the fourth parameter (the ratio of the X-ray fluence emitted in the prompt 
phase to the fluence from the early afterglow) we use predictions of the standard fireball 
model. The standard fireball model (Wijers, Rees, & Meszaros 1997; Sari, Piran, & Narayan 
1998; Sari & Piran 1999) predicts that in the afterglow phase, the X-ray flux decreases as a 
power law, as the GRB-induced blast waves sweep up surrounding matter and decelerate. 
The index of this power law may change its value on time scales from minutes to hours. In 
this "external shock" model, both the tail of the prompt X-ray emission and the afterglow 
are produced by the blast wave (Wijers, Rees, & Meszaros 1997). The external shock model 
predicts that the amplitude of X-ray prompt emission during the flnal part of the GRB 
should be equal to the back-extrapolated flux from the afterglow. This back-extrapolation 
holds for most GRBs that were observed by BeppoSAX (Frontera et al. 2000a). 

However, because the afterglow measurements by BeppoSAX start 6-8 hours after the 
GRB, there is no direct measurement of the X-ray flux immediately after the GRB. If the 
X-ray afterglow decays as a power law, its integral fluence can reach 50%-200% of the prompt 
emission (Frontera et al. 2000a). We have assembled from the literature decay indices for 15 
GRB, and have calculated indices for 10 more, based on the published estimates of prompt 
and afterglow fluxes. The decay indices a fall in the range -1.0 to -1.9, most commonly 
about ~1.3, and the distribution of indices can be satisfactorily fit with a Gaussian. 

In Fig. 5, we combine all four parameters discussed above to obtain a log N-log S dis- 
tribution for X-ray emission from GRBs. This estimate is based on the gamma-ray (50-300 
keV) log N-log S distribution for all BATSE GRBs with duration T90 > 2 s. The solid curve 
is based on our XDGRB distribution of Rx/-y , sampled by Monte Carlo methods. To the 
prompt X-ray fluences we add a contribution from the afterglow obtained by Monte Carlo 
sampling of the afterglow decay parameter a and integrating the afterglow out to 10^ sec- 
onds. These curves plot our estimate of the X-ray fluence from cataloged, triggered BATSE 
events. 
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We also show for comparison the expected FXT distributions for less-preferred distribu- 
tions of Rx/y'- the mean and logarithmic mean of Ginga (0.07 and 0.24, Strohmayer et al. 1998) 
and our 'BATSE" distribution. The less-preferred Rxf-y yield numbers of FXTs from GRBs 
that fall below the preferred estimate. 

We would like to understand how these curves would change if we included all the GRBs 
that are too faint to trigger BATSE. Some of them may be, in fact. X-ray rich GRB or XRFs 
(see Kippen et al. 2002). Wc can get some notion of this correction by adding to our sample 
the sample of untriggcred BATSE events obtained by Kommers et al. (1998) and Stern et 
al. (2000), which enhance the BATSE sample of the faint end. Wc find little change for 
all but the faintest part of the XRF log N-log S curve. Fig. 6 shows the XRF log N-log S 
distribution for triggered, untriggered, and all known BATSE GRBs, based on the XDGRB 
estimate of Rxf-y- 

Fig. 7 compares our prediction of the FXT's from GRBs to our best-fit distribution of 
all FXT. The plotted distribution (solid curve) is based on the XDGRB estimate of Rx/j, 
and the population of all (triggered plus untriggered) BATSE GRBs. 

The distribution derived from the 'BATSE'i?x/7 (shown with dotted curve). BeppoSAX 
data (Heise, in't Zand, & Kuulkers 2000, Heise et al. 2001) suggest that about 70% of FXT 
detected by BeppoSAX are GRB-related. The 'BATSE' Rx/^ would yeld a value of 0.1-1% 
for this ratio, far less than observed. 

Fig. 7 also compares our estimate of the log N-Iog S curve for FXT from GRBs with 
the rate of candidate GRB afterglows reported by Greiner et al. (2000). For any reasonable 
assumption about the Rxf-y distribution, the Greiner et al. (2000) rate of events without 
optical counterparts is in good agreement with our prediction, suggesting that a substantial 
fraction of their residual, non-stellar, events are indeed GRBs. 

It is clear that the fraction of fast X-ray transients that can be attributed to classical 
gamma-ray bursts drops with decreasing fluence. The predicted X-ray event rate for "nor- 
mal" GRBs falls a factor of ~ 3-5 below the total FXT rate at high fluences, but a factor of 
100 below at the faint end. Even if we are off by a factor of two or three due to unknown 

effects, a large majority of FXTs are not counterparts of classical GRBs (e.g, not X-ray 
flashes), especially at the faint end. However, we cannot rule out a possibility that exists a 
population of X-ray rich GRB-like events that comprise a large fraction of the FXTs at the 
faint end. 
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5. OTHER SOURCES of FXTs 

5.1. Possible candidates 

The timescales and spectra of FXTs are diverse. Timescales range from seconds to hours, 
and spectra range from very hard (kT>20 keV, Rappaport et al. 1976) to soft (kT~l keV, 
Swank et al. 1978). PM83 suggested the time profiles could be divided into classes, perhaps 
indicating multiple progenitor classes. The few FXT identifications made to date confirm 
the heterogeneous nature of FXTs. 

The lower right corner of our X-ray log N-log S curve is formed by several events with 
duration up to one day, which were seen mostly by WATCH (CT99). A T Tauri star was 
suggested by CT99 as the progenitor for one of the events, while the rest were probably 
generated by X-ray binaries. Several X-ray binaries have been recently found to gener- 
ate short, powerful outbursts that might be classified as FXTs. For example, outbursts 
from V4641 Sgr up to 12.2 Crab (2-12 keV) with durations less than a day were observed 
by the RXTE and BeppoSAX (Smith, Levine, & Morgan 1999, (in't Zand et al. 2000b; 
Wijnands & van der Klis 2000)).. During the RXTE PCA observation, V4641 Sgr showed 
X-ray fluctuations by a factor of 4 on timescales of seconds and ^ 500 on timescales of 
minutes (Wijnands Sz van der Klis 2000) . The mass function identifies the system as a high- 
mass X-ray binary harboring a black hole (Orosz et al. 2000). A slightly longer type of X-ray 
transient, CI Cam, with a single short outburst about 2 days long, was also discovered by 
RXTE (Smith et al. 1998, Revnivtsev, Emel'yanov & Borozdin 1999). 

"Super type I X-ray bursts" that last from 30 min to 3 hours were reported recently 
from several low-mass X-ray binaries (e.g., Cornelisse et al. 2000; Wijnands 2001). Similar 
events might been detected as high-fiuence FXTs by low-resolution experiments. These or 
similar sources might contribute to the lower right corner of the FXT log N-log S curve. 

In't Zand et al. (1998), and Kaptein et al. (2000) announced the discovery of type 
I X-ray bursts from 'empty' places. These bursts, observed from locations without steady 
emission, would be naturally classified as FXTs. The origin of such bursts might be low 
accretion rate neutron star binaries. We can get an idea of the rate of 'empty' place bursts 
from BeppoSAX(in't Zand 2001). It observed the Galactic Center region for Tq^s = 4 x 10® s 
and detected about ~1500 X-ray type I bursts from A^6=31 sources. At the same time 
it observed Be— 4: X-ray type I bursts from 4 'empty place' X-ray bursters. If all these 
events were generated on low-accretion X-ray neutron stars (with, say. X-ray luminosity 
less than several percent of the Eddington hmit), and their space distribution follows the 
space distribution of normal X-ray bursters, we can estimate the total number of bursts 
from 'empty place' per year per sky as: Ngky — (iyear/Tobs) ■ {Nt/Nf,) ■ ~ 50 events, where 
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Nt ~ 50 is the total number of known normal bursters. If their distribution is isotropic, the 
number of bursts from 'empty' places would be at least an order of magnitude greater. This 
simplified estimate agrees well with the BeppoSAX estimate (Cornelisse et al. 2002a). They 
found that the total number of such sources in our Galaxy is between 30 and 4000, while the 
recurrence time between the bursts is a half year (for 30 sources) and several tens of years 
(for 4000 sources) The fiuence from such events would be comparable with outbursts from 
RS CVn binaries, yielding about 50-100 events per sky in the 10 — 1000Cra6 x sec region. 

PM83 identified a large fraction (6 out of 27) of Ariel-5's fast X-ray transients with 
RS CVn binary systems. Another RS CVn object was proposed as a possible FXT source 
by AW86. RS CVn systems are binaries formed by a cool giant or subgiant with an active 
corona and a less massive companion in a close synchronous orbit. RS CVn binaries typically 
exhibit a peak flux ~10^^ erg/s for a duration of 1 to 10 h. Outbursts from nearby RS CVn 
systems would have a fiuence that puts them in the middle part of our diagram. This is 
certainly true for the identifications proposed by PM83 and AW86. 

Our analysis of GRB X-rays indicates that this contribution to the FXT log N-log S is 
most significant around 100-1000 Crab x sec. Heise et al. (2000) reported that about 70% of 
all FXTs detected with BeppoSAX are GRB-related, i.e., either conventional GRBs or XRFs. 

1 Our best estimate (XDGRB) indicates that the GRB contribution to FXTs (sec Rx/'y 
distribution discussion) peaks at 20-30%. However, taking into account the uncertainties 
of both studies, we consider these results as moderately consistent between themselves, 
and consistent with Grindlay 1999. We note that BeppoSAX/WFC detected ^1.5x10^ 
thermonuclear X-ray bursts from ~35 sources, but did not classify them as FXTs because 
they had clear identifications. Some of these bursts might have been identified as FXTs 
by earlier experiments like Ariel-5. This would increase the total number of FXTs and 
decrease the GRB fraction. Our estimate of GRB-related FXTs is also consistent with the 
ROSAT result (Greiner et al. 2000) at low fluences. We have shown (see Fig. 7) that the 
contribution of classical BATSE GRBs drops dramatically for lower fluences. Other sources 
must dominate FXT statistics that region. 

C86 found 6 of their 10 FXTs to be flares from late type dMe-dKe stars, i.e., M or 
K dwarfs with Balmer lines in emission. AW86 also identified 3 of 10 HEAO-1 Al FXTs 
with fiare stars. Flares of dKe-dMe have peak fiuxes of 10^^-10^^ erg s~^ and duration from 
~minutes to ~hours. Their spectra have temperature ~1 keV, and they should be well 
detected by ROSAT. Indeed, V98 found star-like counterparts for 132 out of 141 reported 
FXTs. We therefore suggest that the least luminous and most numerous part of log N-log S 
distribution is formed by nearby fiare stars. 

There is however a fraction of lower-fiuence FXTs that are not identified with star-like 
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objects. While the contribution of classical GRBs should be neghgible in this region, other 
types of progenitors may play a more important role. One of them, gamma-ray quiet GRBs 
might arise in several ways (e.g. MacFayden & Wooslcy 1999; Meszaros & Gruzinov 2000). 
Recent BeppoSAX results have shown that some GRBs arc faint in the standard 50-300 keV 
BATSE band but bright in X-rays (Kippen et al. 2001; Kippen et al. 2002). Out of the 53 
GRB-like transients detected with the BeppoSAX WFC, 17 were not registered by Bep- 
poSAX GRB monitor, which is sensitive in 40-400 keV energy region. Kippen et al. 2001 
found that 9 of such XRFs were also seen by BATSE and recorded as untriggered events at 
BATSE database. According to Sazonov et al. (1998) 10% of WATCH events had detectable 
emission in 8-20 keV energy channel but no significant emission in the hard channel. Our 
analysis (see Fig.6) confirms that untriggered BATSE events contribute only modestly to 
the lower fluence part of log N-log S if they follow the same Rx/^y distribution as the trig- 
gered events. But GRB-like outbursts with softer spectra might be more significant at lower 
fluences. 

Other exotic progenitors may contribute also. Schafer, King, & Deliyannis 2000, sug- 
gested that ordinary main sequence stars like our Sun may exhibit rare powerful superfiares 
with 10'^'^ - 10'^*^ ergs, some of them in X-rays. However, one cannot quantify their contribu- 
tion until one better understand the frequency of such superfiares. 

Extragalactic sources may also generate FXT phenomena. X-ray variability is a fun- 
damental property of Active Galactic Nuclei. In the last decade it was found that BL Lac 
objects (blazars) exhibit strong correlated variability in both X-ray and TeV gamma-ray 
bands on short timescales of days to hours (Catanese et al. 1997; Marachi et al. 1999). The 
brightest of such events may be detected as FXTs. 

5.2. Contributions to the Fluence Distribution 

Previous searches for FXT progenitors demonstrated that we probably deal with a 
mix of close (flare stars) and extragalactic (GRB-related) events. The sky distribution 
for both these populations is expected to be isotropic. This is consistent with the sky 
distribution of detected FXT (Fig. 3). However, the absense of reliable exposure maps for 
many observations and possible selection effects leave room for the existence among FXTs 
of a non-isotropic population with the Galactic disk or bulge distribution. The slope of the 
log N-log S distribution of such additional component is expected to be flatter than -1.5 for 
a variety of possible progenitors. 

The most natural candidates for an additional component of FXTs are low mass X-ray 
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binaries (LMXBs). Cornelisse et al. 2002b have analyzed BeppoSAX data and found several 
X-ray type-I bursts from 'empty' places, where no persistent emission was detected down 
to a several mCrab limit. Cocchi et al. 2001 have proposed a new class of low- luminosity 
bursters. Cornelisse et al. 2002b have shown that all detected bursts from 'empty' places 
are indeed the members of such class and follow LMXB space distribution. Since LMXBs 
are concentrated in the Galactic Bulge (see e.g. recent review by Grimm, Gilfanov, & 
Sunyaev 2002), any associated extra FXTs should be clustered around the Galactic Center. 
However, most early surveys had difficulty in detecting FXTs from this part of the sky 
(Warwick et al. 1981). Their localization accuracy was poor (typically several degrees), and 
any X-ray flash from a populated region of the sky would tend to be attributed to known X- 
ray sources. Hence, the population of LMXB with low persistent flux and emitting rare bursts 
would need to be be quite numerous to have been clearly detected by those instruments. 
Emel'yanov et al. (2001) have shown that TTM observations of the Galactic Bulge put strong 
constraints on the number-frequency function of low-accreting X-ray bursters indicating that 
the number of low-luminosity X-ray bursters x the burst frequency from a burster is small. 
If the FXTs that are distributed as a bulge population are fewer than several percent of all 
FXTs, they cannot be found with the usual statistical tests (e.g., a K-S test would not identify 
a bulge fraction less than 3% of the total as signiflcant) . The deviation from uniformity in 
the FXT distribution would be further masked by the poor ability of previous studies to flnd 
FXTs at low latitudes (see e.g. Warwick et al. 1981, for the discussion of source confusion 
in Ariel-5 survey). 

Our two most reliably identifled source populations, gamma-ray bursts and stellar flares, 
although isotropically distributed, might yield a log N-log S distribution with a slope signifl- 
cantly flatter than -3/2 (as would be expected for a uniform distribution of standard candles). 
In the case of gamma-ray bursts, the distribution is known to flatten due to cosmological 
effects. Even more important than the cosmological flattening, the luminosity distribution 
of FXTs from GRBs is flattened by the broad distribution of Rx/R-y- We demonstrate below 
that nearby flare stars would also produce a distribution flatter than -3/2. 

EXOSAT observed 22 late-type flare stars (dMe-dKc) in the vicinity of the Sun (2-20 
pc) in the 1-10 keV energy region. They found that these stars generated X-ray flares with 
total energies from 10^° to 10^^ ergs and A'"(> E) ~ E~^ '^ (Pallavicini, Taghaferri, & Stella 
1990), e.g. with a differential frequency-flare distribution fits a power law N(E)~ E"^-"^ . A 
similar distribution function has been observed for the Sun, with differential frequency-flare 
energy distribution index F = 1.8 (Gershberg 1989) or 2.0 (Veronig et al. 2002), and for 
optical flares on the M dwarf flare star AD Leo (Pettersen, Coleman & Evans 1984). 

The integral distribution of X-ray flares from flare stars is: 
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where n(r, E) oc r"^ x E^"' is the differential distribution function of X-ray ffares from 
ffare stars, E is the ffare energy, Rmin is the nearest ffare star, and S=E/(4 vrr^) is the 
measured fluence. This distribution holds for sources uniformly distributed through the 
space (as is certainly within a radius of 100 pc), with a differential frequency- flare energy 
function described with power law E oc E~'^ . We derive: 



oo 

j E-^dE - {47rr^S)-^+^ oc r^-^^ x S^-^ (3) 



B=47rr25 



substituting 7 = a + 1, where a is the slope of the cumulative frequency-flare energy 
function, we get: 



00 

N{> S)(x X J r^-^'^dr (4) 



If the total flare energy dynamical range E — E^''[Efnin, E^ax], is greater than the 
distance dynamical range R^[R — {Rmin, R-max)], then for a range of detected fluences S — 
[Emin/ Rmim ^max/ Rmaxl slope of the detected Cumulative number-fluence function (A^(> 
S)) will be the same as the slope of the integral frequency- flare energy function A'"(> E) oc 

EXOSAT data confirm that A^(> E) oc E'^-'^ for at least Emax/Emin = 10"^ (Pallavicini, 
Taghaferri, & Stella 1990). So, for example, if Rl/R^ — 100, the dynamical range where 
A^(> S) oc S-^-^ is S1/S2 = 100. 

Integrating over the full range of distances (2-100 pc) and flare energies (10^^-10^^), we 
estimate that the integral fluence distribution (A^(> S)) from stellar flares has a slope that 
steepens slowly from —a to -3/2 with increasing fluence (Fig. 8). Note that over 3 orders of 
magnitude in fluence, it is close to a power-law with slope ~ —0.6 ~ —a. 

The intrinsic distribution of flare energy drives the observed log N-log S fluence distri- 
bution for the stellar component of FXTs. The fluence (Crabxs) in the 2-10 keV region is 
related to the flare energy as 3 x 10^° (erg) x M (at Crab x sec) x D'^ (at pc^). The HEAO-1 
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experiments had a fluence detection threshold about 1 Crabxs, while Ariel- 5 had a threshold 
of about 100 Crabxs. As a result, HEAO-1 detected flares with energy L > 1 x 10^^ ergs 
from sources located within 2 pc, and the most energetic flares from distances up to 100 pc. 
Ariel-5 was able to detect flares with L > 1 x 10^^ within 2 pc, and assuming the most ener- 
getic flares within 10 pc. Taking the star density near the Sun ~0.1 star pc"^, and assuming, 
based on the standard 4-component model of the Galaxy (Bahcall & Soneira 1986), that stars 
are distributed uniformly within 100 pc of the Sun, HEAO-1 should have seen flares from up 
to 4 ■ 10^ X Kfs flare stars, while Ariel-5 was able to detect 400 x Kfg flare stars. The catalog 
of chromospherically active stars (Strassmeier et al. 1993) includes 3 stars closer than 10 pc, 
while Pallavicini, Tagliaferri & Stella (1990) observed 12 flare stars within 10 pc. Basing on 
these data we accept Kfs = 2.5 • 10"^ as a lower limit for fraction of dMe-dKe stars around 
the Sun. In this case the HEAO-1 experiments were able to detect X-ray flares from 10"' flare 
stars. Because the HEAO-1 experiments were sensitive to fainter events than Ariel-5, and 
because flare star FXTs arc generally fainter than RS CVn FXTs, we expect that the ratio 
between late-type flare stars and RS CVns should increase as we move from the Ariel-5 to 
the HEAO-1 fluence band. This is consistent with their proposed identiflcations: HEAO-1 
identifled 10 events with dMe-dKe stars, while Ariel-5 reported only one such identiflcation. 

In summary, over a wide range of fluences, the intrinsic distribution of flare energies 
dominates over the geometric term, so that the -3/2 power law expected for standard candles 
does not apply. The shallow slope of the EXT distribution is, therefore, obtained naturally 
even for uniformly distributed progenitors. 

6. SUMMARY 

We have produced a composite log N-log S relationship for fast X-ray transients, by 
combining data from several experiments. This relation spans 4 orders of magnitude, and 
is flt to a power law with a slope of — l.O^Qig. Extrapolation of this slope to lower fluences 
is consistent with the sensitive ROSAT and Einstein experiments. The sources of fast X- 
ray transients are undoubtedly heterogeneous. The a ~ -1 power law must come from the 
summation of contributions from several different progenitor classes. Both GRBs and flares 
stars will contribute components to this distribution that arc flatter than the 'isotropic' slope 
of q; = —3/2. The exact shape of the global fluence-frequency relationship might deviate 
from this simple power law, even inside the range wc measured. However, due to the small 
number of events that we analyzed, and resulting large uncertainty in the deduced rate, a 
more complicated luminosity function is not justifled (see for example detailed discussion at 
Murdoch, Crawford, & Jauncey (1973)). 
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We use archival measurements of the X-ray and gamma-ray fluences for 134 GRB to 
estimate the distribution of the X-ray to gamma-ray ratio. This Rx/^y distribution shows the 
presence of highly significant tail of X-ray rich events. Extrapolating from the BATSE GRB 
catalog, we find that the fraction of fast x-ray transients attributable to classical gamma-ray 
bursts is a function of fiuence, being greatest for high fluences and dropping dramatically 
at the faint end. The exact fraction of GRB-associated fast X-ray transients is sensitive to 
our assumptions about the Rxj-y distribution, but our analysis certainly shows that most of 
fainter FXTs are not counterparts of classical gamma-ray bursts. The fraction of FXTs from 
non-GRB sources, such as magnetic flares on nearby stars, is the highest for the faintest 
flashes. 

Our understanding of the fast X-ray transient phenomenon is still poor. The situa- 
tion is similar to that for gamma-ray bursts before BATSE. There have been many FXTs 
detected, but by a variety of instruments each with speciflc restrictions. The data show 
convincingly that FXTs are ubiquitous. But the number of events detected by each par- 
ticular experiment is small, and their statistical analysis is complicated. Wc can identify 
some classes of progenitor sources, but we can determine the fractional contribution of these 
classes only roughly. There remains room for the discovery of new types of sources that 
produce fast X-ray transients, but we certainly cannot yet confirm the existence of new 
classes like "gamma-ray quiet GRBs" or "superflares" at normal stars. A large uniform 
dataset from a dedicated, sensitive wide-fleld X-ray experiment (see e.g. Eraser et al. 2002; 
Priedhorsky et al. 2000; Borozdin et al. 1999) could revolutionize our understanding of the 
EXT phenomenon. The greatest opportunity for discovery may lie in the faintest EXTs, 
posing a challenge to instrument builders. 

The GRB data were taken from the BATSE archive ^. The work of VA was partially 
supported by Russian Foundation for Basic Research grant 01-02-17295. KB is glad to 
acknowledge helpful discussions of the subject with Jean in't Zand and John Heise. VA 
would like to aknowledge the advices from Eugene Churazov and Marat Gilfanov. 

A. Appendix. Selection effects on Rx/'y ratio 

Differences in the sensitivity band and observational strategy signiflcantly affect the 
i?x/7 ratio measured with different instruments. In our analysis we made use of all data 
we aware of to obtain the most instrument-independent estimate possible. We calculated 



^http://gamniaray.nisfc. nasa.gov/batse/grb/catalog/ 
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the mean, median and variance of the measured distributions of Rxj-y for the 'BATSE' and 
XDGRB samples (Table 1) We also fitted the distribution Rx/^ for the XDGRB and 'BATSE' 
samples with a Gaussian, and with a Gaussian plus quadratic polynomial. These functions 
were choosen simply to match the data with a smooth curve. A Gaussian alone cannot fit 
the tail of distribution. The difference between the Gaussian cores of the distribution is in 
the direction expected from the triggering effects, but the tail in the XDGRB distribution 
cannot be caused by triggering effects. The inset on Fig.4 shows the best-fit curves for the 
Gaussian plus quadratic polynomial fit to the both samples, to show qualitatively how they 
differ. 

To understand the effect of triggering on the Rxi-y distribution, we consider the simplified 
case in which the X-ray/gamma-ray fluence ratio follows a log-normal distribution. In other 
words, if a is defined as Fx=exp{a)F^, a is distributed as a Gaussian. The Gaussian 
distribution of a is dN/da oc exp[—{a — /i)^/2(7^], where // is the mean of the distribution 
and (7 is its standard deviation, for events triggered in the band of F^. 

If the log N-log S distribution goes as N S~^, then the X-ray triggered distribution 
of a, for the same population of events, will follow a Gaussian distribution of the same 
width, but centered on ^{Xj-ay) — /J'i'j) + P ■ cr. We assume that the X-ray trigger band is 
the same band as for Fx- In the case of an isotropic, (3 = 3/2 power law, the distribution 
will be displaced 1.5 standard deviations to the right when we trigger on X-rays instead of 
gamma rays. For trigger bands other than those used for Rx or R^, the displacement will 
be different. 

This selection effect may explain the fact that WATCH generally selects softer GRBs 
than BATSE. While BATSE triggers in the 50-300 kcV band, WATCH triggers in the 8- 
100 kcV band. However, the nominal trigger range for Ginga, at 50-400 keV is harder 
even than BATSE. The softer Rx/y ratio obtained by Ginga may be the result of a sec- 
ondary selection done by Strohmayer et al. (1998). While Ginga observed about 120 GRBs 
(Ogasaka et al. 1991; Fenimore et al. 1993), only the GRBs with good quahty spectral data 
in both 7- and X-ray detectors were selected for spectral analysis, possibly driving the se- 
lection of softer. X-ray bright, GRBs. 

We can make a rough check of the selection effect by taking the power law P from BATSE 
measurements, calculating the displacement /3 • cr, and comparing it with the displacement 
between the BATSE and the other, mostly X-ray triggered, Rx/j distributions. We predict a 
displacement for /3 ~ 1.3 (derived from middle part of BATSE GRBs log N-log S) of 0.6-0.7, 

compared to a displacement in range of 0.5-0.7 (Table 1), between the Rx/-y distributions for 
'BATSE' and XDGRB. The direction is as expected: triggering at softer energies moves the 
Rxf-y distribution to larger values. 
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We have tested whether the difference between the 'BATSE' and XDGRB distributions 

was caused by the effect of trigger-band selection on a single distribution of GRBs, or if 
there existed an additional X-ray rich GRB population not detected by BATSE. To do this 
we have fit the -Rx/7 distribution of the BATSE GRB sample to a log-normal distribution. 
A maximum-likelihood estimate of the fit parameters gives fi{Xray) ~ and a = 0.5. This 
result shows that except of the triggering band effect there is also an additional contribution 
of X-ray rich events for XDGRB case ^. 



See also discussion at Kippen et al. 2002 
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Fig. 1. — Maximum-likelihood confidence limits for the parameters of a A^(> S) = AqX S 
power-law distribution of FXTs. The contour lines show 1, 2, and 3 a significance levels. 
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Fig. 2. — Observed log N-log S distribution for FXTs, assuming a power law fluence distri- 
bution {N{> S) = Aq X S~°'). The thick solid line is best-fit power law derived from a max- 
imum likeUhood fit to Ariel-5, HEAO-1, and WATCH data, with a = 0.95 and = 4 x 10^ 
The shadowed areas show the measurements of the log N-log S distribution from each ex- 
periment, including both event statistics and systematic errors. The cross denotes ROSAT 
events (V98), and the diamond Einstein FXTs The dashed and dotted curves show the range 
of distribution that follow from the HEO-1 A2 events, given extreme possibihties for event 
duration. (GHH96) 
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Fig. 3. — Sky distribution of the FXTs included in the log N-log S distribution of Fig.2. 
The crosses are the Ariel-5, HEAO-1, and WATCH events; triangles are the ROSAT and 
Einstein events 
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Fig. 4:.— Rx/^ distribution for GRBs. (2-10 keV fluences /50-300 keV fiuences). The solid 
line is from Ginga, RXTE/Konus, BeppoSAX, and WATCH events (XDGRB case), and the 
dashed hne from BATSE events. The former were measured in the X-rays, but the latter 
were extrapolated from the BATSE band (above 20 keV) to the X-rays. The inset shows the 
best fit to these distributions, using a Gaussian plus a quadratic polynomial (fit parameters 
arc given at Table 1). The solid curve is the XDGRB sample, and the dashed curve is the 
'BATSE' sample. 
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XRF fluence (10 ^ erg cm ^) 

Fig. 5. — The expected log N-log S X-ray distribution (summing prompt and afterglow 
emission) for BATSE triggered GRBs, based on different distributions of the i?x/7 ratio. The 
thick sohd curve is based on the XDGRB distribution Rx/'y, derived from Ginga, BeppoSAX, 
RXTE/Konus and WATCH events (sohd line in Fig.4). The shadowed area shows la errors 
for this solid curve, based on Monte Carlo sampling of the Rxf-y distribution. The dashed 
curve is based on a i?x/7 derived from BATSE spectra ('BATSE' case). The error region 
for the 'BATSE' case has a similar size as for the XDGRB case. Finally, the dotted curves 
shows the expected log N-log S distribution for the fixed values of 0.07 and 0.24 found by 
Ginga (their logarithmic and arithmetic averages; Strohmayer et al. 1998). 
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Fig. 6. — The dashed curve is the expected XRF log N-log S for triggered BATSE GRB, 
with Tgo > 2 sec, using the XDGRB estimate of i?x/7 (as for the sohd hue of Fig. 4). 
The dotted curve is the expected contribution from untriggered BATSE GRB, assuming 
the same Rxj-^ distribution, and the sohd curve is sum of the triggered and untriggered 
contributions. The small difference suggests that undetected GRBs are unlikely to make a 
significant contribution to the FXT distribution. 
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Fig. 7. — Contribution of gamma-ray bursts to the model log N-log S distribution for X-ray 
flashes. The XRF distribution parameters are taken from Fig.l. The thin solid curve is 
the GRB-derived distribution from Fig.6 (XDGRB case), while the thin dotted curve is the 
expected distribution from the 'BATSE' case. The cross is derived from the ROSAT survey 
of afterglow-hke events that were not attributed to late-type stars (Greiner et al. 2000). 
Despite the uncertainties in our analysis, it is clear the classical gamma-ray bursts are more 
important contributor to bright XRFs than to faint XRFs. 
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Fig. 8. — The expected log N-log S FXT distribution for late-type flare stars in the vicinity 
of the Sun. The normalization is arbitrary. The slope for lower fluences is dominated by the 
intrinsic number-energy distribution of stellar flares, rather than the 3/2 slope of standard 
candles in an isotropic geometry. 
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Tablc 1: log{Rx/^) distribution parameters for different models 



Case 


Distribution 






Model parameters 




parameters 




Gaussian 


Gaussian -|- Quadratic polynomial 




Mean Median 


a 


/i a 


/i a A" C 


BATSE" 


-1.67 -1.68 


0.48 


-1.53 0.43 


-1.53 0.43 0.0 0.0 0.0 


XDGRB^ 


-0.95 -1.12 


0.68 


-1.06 0.30 


-1.08 0.22 -0.01 -0.01 0.04 



« BATSE events 
Ginga + WATCH + RXTE + BcppoSAX events, XDGRB case 
A, B & C represent parameters of the quadratic polynomial Ax^ + Bx -\- C 



